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Abstract. The technique of transformation optics (TO) is an elegant method for 
the design of electromagnetic media with tailored optical properties. In this paper, we 
focus on the formal structure of TO theory. By using a complete covariant formalism, 
we present a general transformation law that holds for arbitrary materials including 
bianisotropic, magneto-optical, nonlinear and moving media. Due to the principle 
of general covariance, the formalism is applicable to arbitrary space-time coordinate 
transformations and automatically accounts for magneto-electric coupling terms. The 
formalism is demonstrated for the calculation of the second harmonic generation in a 
twisted TO concentrator. 



1. Introduction 

The field of transformation optics (TO) has drawn a lot of scientific interest in the 
last few years [1-6]. By this design methodology, the form-invariance of Maxwell's 
equations under coordinate transformation is used to tailor the optical properties of 
an electrodynamic medium. The majorities of TO applications focus on the design of 
the linear material parameters. In that coordinate transformation is used to 

engineer the linear constitutive parameters of a medium such that the wave trajectory 
follows a desired path [7]. An interesting application of this concept is the realization 
of electromagnetic invisibility cloaks [8-12] — devices in which light is guided around 
a certain region of space rendering the interior of the region invisible for an external 
observer. Cloaking devices belong to the most prominent applications in which the TO 
concept is successfully applied for designing the linear material properties of a medium 
and have been extensively reported in the literature [13-19]. 

In contrast, only little work in the research of TO media addresses the 
transformation of nonlinear material properties. Media with nonlinear response, 
however, provide a number of interesting effects including sum- and difference- 
frequency generation, parametric amplification and oscillation, stimulated scattering, 
self-phase modulation and self-focusing [20-28] and play a key role in modern optical 
technology [29-31]. Consequently, the extension of the TO concept to nonlinear media 
is expected to offer a variety of new opportunities for engineering optical media and for 
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the construction of novel electromagnetic devices [32]. Apart from that, the TO concept 
also provides a promising calculation method in modeling complex nonlinear systems 
whenever it is possible to find a coordinate transformation such that the geometry of 
the system takes a much simpler form in the transformed space. 

A basic prerequisite for a successful integration of nonlinear effects into the TO 
concept is a general transformation law for linear and nonlinear material parameters 
under space-time coordinate transformation. A first step in this direction is suggested 
in [32] where the transformation of certain classes of nonlinear materials under purely 
spatial transformations is studied. The aim of the following paper is to generalize 
this approach and to provide a rigorous theoretical framework for arbitrary nonlinear 
materials and for both temporal and spatial coordinate transformations. The formalism 
is presented in a manifestly covariant form that allows the simultaneous treatment of 
electric, magnetic and magneto-electric cross-coupling effects and encompasses all types 
of space-time transformations with a particular reference to moving media [33-35]. 

The paper is organized as follows: In the first part, we introduce the used 
tensor notation and explain how the defined tensors are related to the common 
linear and nonlinear material parameters such as the permittivity, permeability or the 
quadratic electro-optic coefficients. Subsequently, we exploit the fundamental principle 
of relativity to derive a general transformation law for nonlinear media under continuous 
space-time coordinate transformations. In this context, a special emphasis is given 
to the transformation of nonlinear constitutive parameters in moving media and the 
class of time- independent, spatial transformations which play a central role for the 
design of TO devices. In the final part of the paper, we illustrate and explain the 
derived expressions by means of an explicit calculation example. By deriving the 
second harmonic generation in a nonlinear, twisted field concentrator, we show that an 
appropriate coordinate transformation can provide a significant alleviation in the formal 
treatment of a complex nonlinear problem and, thus, allows a convenient calculation of 
an otherwise sophisticated process. 

2. The covariant material equation 

In order to establish a common basis for the following discussion and to introduce 
the used notation, we start with a short review of the covariant description of the 
electro dynamic theory. In this paper, all quantities are expressed in Gaussian units. In 
this case, the Maxwell equations in matter take the form: 



div B = 0, 



IdB 



= 0, 



(1) 



rot E + - 



c dt 



div D = 47rp, 



rot H - - 



ldD 



47T . 
= —J 



(2) 



c dt 



c 



where E, H, D, B, p, j and c denote the electric and magnetic field, the electric and 
magnetic flux density, the charge, the current density and the speed of light in vacuum, 
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respectively. By using the four-current (j v ) = (cp,j), the antisymmetric field strength 
tensor 



(3) 
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and the antisymmetric displacement tensor 
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(4) 



the Maxwell equations can be covariantly expressed as: 

d*F Vff + dPF^ 



d"F^ = 0, 
c 



(5) 
(6) 



where we use the Einstein summation conventionf and the contravariant derivation 
= {\d t , —V). In the covariant notation, (5) contains the two homogeneous Maxwell 
equations whereas (6) contains the two inhomogeneous Maxwell equations. Note that 
the signs used in the tensor notation above depend on the convention used for the metric 
tensor. In this paper, we use the metric tensor given by rf v = diag(l, —1,-1,-1). The 
relation between the tensors F^ and is given by 



V = F 

11V - 1 [11/ 



[IV 



(7) 



with the polarization-magnetization tensor 



(?V) 



Tensor equation (7) is equivalent to the common material equations: 

D = E + AttP, 
H = B — AnM. 
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(8) 



(9) 



In many cases, the quantitative relation between the polarization-magnetization 
tensor and the field strength tensor F^ u can be expressed in a power series according 



f Throughout this paper, Greek indices run from to 3 while Latin indices run from 1 to 3. 
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to: 

V =\ aK F + Y CTKa/3 F F r + Y aKa ^ 5 F F rF H 

= E* A ^"^A-.^- (10) 

n>l 

Note that in this notation, the index pair \xv is fixed while the other indices are summed 
over. From the antisymmetry of the tensors F^ v and and the commutativity of 
the products F^^F^^, it follows that the coefficients xZlf 1 ' an ^ n are antisymmetric 
under exchange of //■<->■ v, symmetric under exchange of two pairs aij/^ <->■ (Xj/3j and 
antisymmetric under exchange of two indices on -B- within one pair For the 

quadratic term, this means for example: 

ai/3iQ!2/32 _ v ai/3ia2/32 _ v «2/32ai/3i _ v /3iaic»2/32 fill 

A^u Avfi a^iv Afiv ' V / 

It is obvious that additional, inherent symmetries (such as the Kleinman symmetry or 
spatial symmetries given by the point symmetry class of the medium) further reduce the 
number of independent coefficients. However, in order to provide a general description, 
we only consider the basic symmetries given in (11). Inserting (10) into the material 
equation (7) yields: 

V = F + 4ttV 

r iiv \ / J Afiu r aiPi ■ ■ ■ r a n p n 

n>l 

= 47r^ X ^-^F Ql/3l ...F a ^ (12) 

n>l 

where in the last line, we performed a re-definition of the linear coefficient x^lf 1 to 
include the free-space space contribution (see appendix). 



2.1. The linear term 

To become familiar with the tensor notation, it is instructive to first consider only the 
linear term on the right-hand side of (12): 

V^ = 4n X ™F aK . (13) 

By using the symmetry properties of it follows for the components of the electric 
displacement field D: 

A = V 0l = 4vr X ™ F aK = 8ttx2' F 0j + 4vr X ^ F kl 

= 4Ej + $B 3 (14) 

where e\ denotes the permittivity tensor and ^ is the tensor of the electro-magnetic 
coupling (remember, Latin indices run from 1 to 3). The last line follows from the fact 
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that F j addresses the components of the -E-field while F ki addresses the components of 
the £?-field. In a similar way, one finds for the components of the i3"-field: 

TT \n mn T) 9-7rn mn -\/ CTK TP 

= ({e, + (// ■ )'i>j (is) 

with the totally antisymmetric Levi-civita tensor g™ n (with gf 3 = 1), the permeability [i\ 
and the magneto-electric coupling tensor Q. In summary, the linear tensor equation 
V^ u = 4:irxZv F aR can be expressed in the common three-formalism as: 

/ n \ / , f \ I tt<\ 

(16) 

This is the general equation of a linear bianisotropic medium. The relations between 
the parameters e, //, £ and ( and the four-tensor can be derived by comparing the 
coefficients occurring in (14) and (15). As shown in the appendix, one finds: 

4 = snxt el = -^ginnxT, 

CI = -4^r n xL, (^Yi = tonLufxW- (17) 

2.2. Nonlinear terms 

Next, we consider quadratic contributions to the polarization- magnetization tensor V^, 
i.e. contributions that depend on the product of two components of the field strength 
tensor. These are given by the second summand in (10) and have the form: 

T>(2) _ uKafi t? t? 
' ixv /V/ii/ 1 o-k 1 a/3- 

For simplicity, we restrict to the components of the second-order electric polarization 

(2) (2) 

P> = (a similar derivation applies for the magnetization). As for the linear term, 
we can use the symmetry properties of 1° split the summation into three terms: 

p(2) = ™c*P p p 

A-.0k0m 771 rn i /K.Ofcmn 771 771 i ,.fc£mn 771 771 

a£EjEu + + cf 5^ (18) 

Pockels effect, multi-wave mixing Faraday effect 

where we have introduced the nonlinear material coefficients a{ k , ty 1 *, c{ k . In other 
words, the expression vj?J = F aK F a p describes all second-order electric, magnetic 

and magneto-electric cross-coupling effects in a single equation. Similarly, one finds for 
the nonlinear polarization of the third order: 

p (3) = FaKFaf}Fiw 

-• jkl EjE k Ei + \? kl EjE k Bi + ^%BuB^ + dfB^Bx (19) 

Cotton-Mouton effect 
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and so on. Note: while the quadratic contributions in (18) contain 6 + 9 + 6 = 21 
independent coefficients, the cubic contributions in (19) contain already 10 + 18 + 18 + 
10 = 56 independent coefficients (see appendix). 

3. Coordinate transformation 

The particular advantage of the covariant formulation is its form-invariance under 
coordinate transformations 

x a — -> x a '(x a ) (20) 

where (x a ) = (ct,r) is the coordinate vector of the space time. By = dx a ' /dx a 
and \A\ := det(A°,) we denote the Jacobian matrix and its determinant. To clarify the 
following discussion, we introduce the convention that, if a product of A^ occurs in a 
transformation formula, the kernel symbol A is written only once, e.g. 

= (21) 
The electromagnetic fields F M „ and in the primed and unprimed coordinate systems 
are related by [33, Chap. 3.2]: 

V^, = {A^A^V^. (22) 

The Maxwell equations (5) and (6) have the same form in the primed coordinate system 
as in the unprimed system due to their natural form-invariance [33, Chap. 3.2], that is: 

d^F v , & + VF^ + d v 'F clll , = 0, 

= ( 23 ) 

with j v i = \A\~ 1 A U U , j v . In order to achieve form-invariance also for the constitutive 
relation (12), the material coefficients must transform as: 

^[P'l-a'nP'n _ I A\-l A» v A a lPl- a 'nP'n ^.OixP^-a-nPn (0A\ 
A-illv' — \^\ 71 M V ^aiPl-OtnPn X-^U ■ 

The proof of (24) is not difficult. To see this, we multiply both sides of (12) by \A\^ 
insert in the right-hand side the identity (see [33], Chap. 1.4) 

A a 'lP'v a 'nPn AOCl/3l—a n Pn _ AOCl/3l---a n p'n _ 1 

and replace the unprimed quantities with the corresponding primed ones. This yields: 

n>l 

1—1 At* v A a lP'l"' a nP'n 



I 4 1-1 Ai iu Tl — Att \ ^1/11 — 1 Af- V A a lPl'" a nPn ^.aiPl-at n p~n A a lPl- a nPn T? T? 

\ A \ \' V ' U ^ - 47r 2^ A M'^' A ai/5i-«n/3n X/w A a'P[-a' n p'J^ ■ ■ ■ 1 ^Pn 

^ » ' „>1 V v '> „ ' 

F a ,a,...F a ,g, 

^ Y / / 1 p l "nPn 

= 4tt J] xffi'^Fw . . . F a , A . (26) 

n>l 
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Thus, the constitutive equation transforms covariantly as required. Relation (24) 
represents the general transformation law for linear and nonlinear materials parameters 
(including complicated magneto-electric coupling terms) and is equally valid for 
arbitrary space-time coordinate transformations. 



3.1. Moving media 

A typical physical situation where coordinate transformations play a role occurs when 
a medium moves relative to the observer. In such cases, the transformation law allows 
the calculation of the material parameters in the reference frame of the observer if the 
material parameters are known in the rest frame of the medium. Consider, for a example, 
a linear, homogeneous medium that is isotropic in its rest frame with the constitutive 
equations: 



D = eE, 
H = p~ x B. 



(27) 



According to (17), the covariant material tensor x„£ of an isotropic medium is related 
to the material parameters e and \i by the equations: 



f -1 = 87rXH 



(28) 



(modulo interchange of indices, e.g. e = — 8nXio, etc.) while x^u i s zero otherwise. By 
means of the transformation rule of (24), we can now calculate the material equation 
in a reference frame in which the medium moves uniformly, say along the x-direction, 
with constant velocity v. The corresponding coordinate transformation is given by 
= A"!,x a where the Jacobi matrix A"' describes a Lorentz boost in the minus x- 



x 



direction: 
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V 



\A\ 



(29) 



with (5 = v/c and 7 = (1 — /3 2 )" 1 / 2 . The inverse of the Jacobi matrix A", is obtained 
by replacing j3 by — j3 (corresponding to a velocity reversal). The material equation in 
the frame in which the medium is moving is obtained by applying the transformation 
law (24) to the material equation (13) (expressed in primed coordinates), that is 



£W = 4tt X ^ F a , K , 



a Ativan 



Frr'K 



(30) 



After a straightforward summing over the indices and by using the relation between the 
field tensors X^v and F a i K > and the E' ', B', D' and H' fields (see Eqs. (3) and (4)) as 
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well as the relation between the material tensor and e and ji given by (28), this can 
be rewritten in the more intuitive form: 



D' = eE' + 7 2 (e - £) " x [b 1 - V - x E') , 



H' = 



(31) 



in agreement with the findings in [36]. Since these equations are form-invariant under 
spatial rotation of the coordinate system, they are equally valid for arbitrary directions 
of the medium's velocity v. Obviously, a linear medium that is isotropic at rest becomes 
bianisotropic if it is moved relative to the observer [35-38]. 

The mixing of the electric and magnetic response in moving media is a general 
effect which can also be observed for the nonlinear material coefficients. Consider, for 
example, the quadratic term of the material tensor given by 111 the frame in 

which the medium is moving the tensor transforms as 



XpJv' a k a I 



• V v crua/z 
fi'u' Xfxv 



(32) 



As before, we can evaluate the sum, collect corresponding terms and identify the 
new nonlinear coefficients. However, in this generality the final expressions get very 
complicated and unintuitive. Therefore, in order to illustrate only the principle of the 
mixing of higher-order susceptibilities in moving media, we restrict to the second-order 
contribution to the electric polarization P- 2 ^ = Xoi^ F aK F a p. As described in (18), 
the relevant components of the material tensor can be further split in terms with Xo?°\ 
X^ kl and x 1 Qi lkl ( m odulo interchange of indices) to separately account for contributions 
that are proportional to EiEj, EiBj and BiBj, respectively. If we now apply the 
transformation (32) separately to these three coefficients, the result can be expressed in 
a compact matrix form: 



(33) 



Obviously, the time-dependence Aq = c~ x dx % ' ' jdt ^ due to the movement of the 
medium implies nonzero off-diagonal elements which induce a mixing of nonlinear 
material coefficients of electric and magnetic nature. This means, for example, that 
a Pockels medium at rest (xoi° l 0) can display a Faraday effect (Xo>i>' k ' 1 ' 0) if ^ is 
moved relative to the observer and vice versa. 



/ ...O'n'O'Z' \ 






4<r 


A ' ' \ 

A mk 1 




/ A .OnOZ \ 
Xoi 1 


^n'n'k'U 
Xo'i' 


aO i n'V 
— A 0'i'n I 




4<r 


AO'k' 
k 




AOi 


\A0V / 




VK'o k ' 


±Kk 


A m'k' , 
^-m k / 




~.,mnkl , 
\A0i / 



Covariant Description of Nonlinear Transformation Optics 



9 



3.2. Time-independent transformations 

We now focus on the special case of time-independent, spatial transformations, 
i.e. transformations with 

t' = t 

x' = x'(x,y,z) 
y' = y'(x,y,z) 

z' = z'{x,y,z) (34) 

which are of particular interest for the design of TO devices with tailored optical 
properties. In this case, we have 

< = = # 4! = g = «. (35) 

Consequently, if the term Af'T a ip> occurs in an expression, we can replace the four- 
index a' by the three-index i' (which runs from 1 to 3) since Af = 0. For example, the 
permittivity transforms under these conditions as 

4 = Bn Xw 



^\AY l ^A\ t bl4x 



= \AY'A%^. (36) 
Accordingly, the permeability and magneto-electric coupling terms in (17) transform as: 



4 = \a 


~ 1 A ij 'n j 


# = \a 


-lAij'd 




Si • 



(37) 

Note that for purely spatial transformations there is no magneto-electric cross-mixing. 

By the general transformation law (24), we can now calculate the transformation 
behavior of the nonlinear material coefficients under spatial coordinate transformations. 
For instance, we exemplarily calculate the second order susceptibility of the electric 
polarization a{ k = Ax^° k which describes nonlinear optical effects such as second- 
harmonic generation or three- wave-mixing. With the help of (24), the general 
transformation of the second order tensor Xu^ * s 

a'da'P' _ I a i-l a (iv a'da'P' ana/} /qo\ 
X/x'u' — I 71 ! ^n'l/anaP Xjw \ M ) 
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and for the special case of a spatial coordinate transformation (i.e. with the conditions 
given in (35)), the second order susceptibility transforms as 

j'k' A O'j'O'k' 

— 41 4 I- 1 /t^O'j'O'fc' ana/3 

= M A \~ lA i'Hk xT k 

= \M- l A\^a? (39) 

where we applied similar calculation steps as in (36). This shows that as soon as 
the relation between the material coefficient of interest (linear or nonlinear) and the 
general covariant tensor xffi 1 '" an/3n is found, the transformation law immediately follows 
from (24). 



4. Twisted, nonlinear field concentrator 

In the following, we demonstrate that nonlinear problems in complex media with 
sophisticated linear and nonlinear optical properties can take a much simpler form in 
an appropriately transformed space. 

To provide an illustrative example of this calculation method, we consider a 
nonlinear, inhomogeneous material that is illuminated by a strong laser field. We assume 
that the polarization of the fundamental wave and the nonlinearity of the material match 
the condition for second harmonic generation (SHG) and, as a goal, we want to calculate 
the spatial field distribution of the SHG wave inside the medium. If the material is highly 
inhomogeneous, both the fundamental and the SHG wave follow a complicated, distorted 
trajectory through the medium which generally hampers a numerical calculation and 
reliable prediction of the SHG progress. However, if it is possible to find a coordinate 
transformation such that the wave propagates uniformly along straight lines in the new 
coordinate system, the wave equation can be readily solved by an ordinary integration 
along the field lines. A subsequent back-transformation then yields the SHG field in the 
physical space. In the following, we demonstrate and explain this calculation method 
for a specific example. 

As a hypothetic, inhomogeneous material, we consider the special case of a TO 
medium, i.e. an artificial material whose permeability and permittivity tensors were 
obtained by applying a coordinate transformation to a homogeneous, isotropic space. 
We suppose that the optical properties of the homogeneous space are similar to that of 
vacuum with a permittivity and permeability equal (or close) to unity. In the following, 
the coordinates of the uniform, isotropic space are indicated by unprimed indices while 
the coordinates of the physical space of the inhomogeneous material are indicated by 
primed indices. 

For the transformation between the primed and unprimed system, we consider the 
following transformation (expressed in cylinder coordinates using r 2 = x 2 + y 2 and 
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= arctan(y/:r)): 



0' 




2 "^^ife" 



0<r < R 2 
R 2 <r< R* 
otherwise 

< r < R 3 
otherwise. 



(40) 



As illustrated in Fig. 1, the transformation compresses space of a cylindrical region with 
radius R 2 into a region with radius R\ at the cost of an expansion of space between 
Ri and R% [15]. In addition, the space experiences a radius-dependent twist about 
the z-axis [39]. Note that the transformation addresses only the x- and y-coordinates 
while the z-coordinate remains unchanged. For this reason, we can restrict the following 
discussion to the xy-submanifold of the space-time manifold. The corresponding Jacobi 
matrix is: 



dr'' 



Ri 

R 2 



«3< 



R3-R1 
R3-R2 

4.R3 blI ^i?3- 



diag(l, 1) 



< r < R 2 



R 2 <r<R 3 



otherwise. 



(41) 



And the determinant of A g is: 




< r < R 2 

R2<r<R 3 

otherwise. 



(42) 



The transformation (40) is expressed in cylindrical coordinates. For the calculation of 
the SHG, however, it is more advantageous to formulate the transformation in Cartesian 
coordinates in the form x' = x'(x,y) and y' = y'(x,y). This is achieved by applying 
an intermediate transformation r = \J x 2 + y 2 and <p = arctan(|//a;) in the unprimed 
system and an inverse intermediate transformation x' = r' cos <f> and y' = r' sin 0' in the 
primed system. With these expressions, the transformation (40) can be re-expressed in 
Cartesian coordinates by applying the following three subsequent transformations: 

/ 



(x,y) — > (r,0) — y (/, 0') — > (x',y'). 



(43) 



According to the chain rule in higher dimension, the Jacobian matrix of a composite 
function is just the product of the Jacobian matrices of the composed functions, that is 



Afogoh = AfA g A h . 



(44) 
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Figure 1. Visualization of the twisted cylindrical concentrator. The mesh grid 
indicates lines with constant values of x and y plotted in the coordinate system of 
x' and y'. The radius of the three displayed circles are R\ — 3, i?2 = 4.5 and Rs = 6. 



As shown in the appendix, the Jacobian determinants of the intermediate 
transformations are \Af \ = 1/r and \Ah\ — r', respectively. Consequently, the Jacobian 
determinant of the composite is 

\A f ogoh\ = \A f \\A g \\A h \ = —\A g \ (45) 



r 

with the Jacobian determinant \A g \ given in (42). Once the Jacobian matrix and its 
determinant are known, we can easily transform the fields and material parameters from 
one coordinate system to the other. 

We now intend to calculate the SHG wave generated in the concentrator. For 
this purpose, we suppose that the material used for the construction of the concentrator 
exhibits a non- vanishing nonlinear susceptibility. To simplify matters, we further assume 
that this nonlinearity obeys the phase matching condition for frequency doubling if the 
fundamental wave and the second harmonic wave are both polarized in the ^-direction. 
In our notation, the corresponding tensor component of the nonlinear susceptibility 
is a' := 4xq^ ' 3 ' (see (18)). For the spatial distribution of a', we suppose that the 
nonlinearity is only present in the inner cylindrical region of the material according to: 

/, /\ U 0<r'<f?! 
a'(r') = <^ (46) 
I otherwise. 

This could, for example, be realized by doping the center of the concentrator with some 
nonlinear material. 
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As fundamental wave we assume an incident monochromatic plane wave that 
initially propagates along the x-axis while the electric field vector is polarized in the 
^-direction. In the physical space of the inhomogeneous material (spanned by primed 
coordinates), the fundamental wave takes the form: 

KW,V',t)=£Le*<>** + W- ut) (47) 

where E'^ denotes the amplitude of the wave and k' w = (k' ux i, k' wyl ) is the wave vector in 
the medium. For the second harmonic wave, we represent the electric field by 

E' 2w {x', y', t) = ejf, y') e^w^w^ (48) 

with the SHG wave vector k 2uJ = 2k' u (phase matched case). In the slowly varying 
amplitude approximation for the second harmonic wave and the undepleted pump 
approximation (i.e. E' w is constant), the wave equation for the SHG amplitude £' 2w {x! ', y') 
is given by [31]: 

with k = —16Triu 2 /c 2 . This is a partial differential equation in two dimensions 
for the second harmonic field amplitude £ 2u (x', u') which is very difficult to solve in 
general. However, the complexity can be significantly reduced if we apply a coordinate 
transformation to the uniform space (spanned by unprimed coordinates). 

According to (22) and the Jacobian matrix of the transformation given by (44) 
(with suppressed subscript / o g o h), the ^-component of the electric field transforms as 

E' z > — F 's' = A^yF^ = F 03 = E z (50) 

since t' — t and z' — z implies Aq,^, = 6^5^. Consequently, the electric fields of both the 
fundamental and SHG wave remain unchanged under coordinate transformation (but 
the functional dependence on the underlying coordinate system may be different). Since 
the unprimed system is homogeneous and isotropic with a refractive index equal to one, 
the electric field of the fundamental wave in the unprimed system has the form of a 
uniform plane wave propagating in the x-direction: 

E u {x,t)=£ u S k °" x - ut '> (51) 

with constant amplitude £ w = S'^ and k ux = uj/c. From this expression, the electric field 
in the primed coordinate system is readily obtained by expressing x by x = x(x',y'). 
The real part of the fundamental wave in the primed and unprimed system is plotted 
in Figs. 2(a) and 2(b), respectively. 

Since the fundamental wave propagates straightly along the x-direction in the 
unprimed system and since we assumed phase matching for the SHG process, it follows 
that the wave vector of the SHG wave is given by k 2uJX = 2k UJX = 2ui/c and k 2u]y = 
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Figure 2. Second harmonic wave generation in a nonlinear cylindrical concentrator. 

(a) Electric field of the fundamental wave in the uniform space (unprimcd system) and 

(b) in the inhomogeneous, physical space (primed system), (c) Electric field of the 
second harmonic wave in the unprimed system and (d) in the primed system. The 
radius of the three displayed circles are R\ = 3, R2 = 4.5 and R3 = 6. 



(on the considered length scale, beam divergence due to diffraction can be neglected). 
Consequently, the wave equation (49) for the SHG process reduces in the unprimed 
system to the simple expression 



-^S 2u} (x,y) = —£*a(x,y) 



(52) 



which can be immediately integrated to: 

X 

hvC 1 

^{x,y) = —£l I dsa{s,y). 



(53) 



The remaining unknown is the nonlinearity a = 4xq3° 3 , i.e. the relevant tensor 
component for the SHG process in the unprimed system. According to the 
transformation law derived in (39), the nonlinearity transforms as: 



/ a.,0'3'0'3' _ A \ a\-1 Ai3'3\0j0k 
— A-3'jkXQi 



MAl-'xT = \A\-'a 



-1 



0'3' 



(54) 
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since z' — z implies that A % y'-^ = S^S^S^. With the determinant \A\ given by (45) and 
(42), the nonlinearity a' defined in (46) and the relation between r and r' given by (40), 
we can calculate the nonlinearity in the unprimed system to be: 

r I otherwise. 

As expected, the nonlinearity in the unprimed system is reduced by a factor of (R1/R2) 2 
compared to the nonlinearity in the primed system (see (46)) because the cylindrical 
region in which the nonlinear substance is located experiences a space expansion from 
radius R\ to Ri if we perform a coordinate transformation from the primed to the 
unprimed system. 

Now we can (even analytically) evaluate the integral (53) for the amplitude S2 U) {x ) y) 
of the SHG wave. A subsequent multiplication with the propagation phasor yields the 
SHG wave in the unprimed system in the form: 

E 2u (x,y,t)=£ 2uJ (x : y)e 2i "W c - t \ (56) 

By applying the inverse transformation (i.e. expressing x and y by x = x(x',y') and 
y = y(x',y')), we finally obtain the SHG wave in the physical space of the primed 
coordinates. The real part of the resulting SHG field distribution in the two coordinate 
systems is plotted in Figs. 2(c) and 2(d), respectively. 

The proposed twisted nonlinear concentrator is certainly a somewhat constructed 
example since the exploited uniform space was already presumed in the design of 
the concentrator. However, the decisive step in the calculation — the straightening 
of the wave trajectories inside the medium by applying an appropriate coordinate 
transformation — is in principle always possible in any inhomogeneous media with 
continuously varying material properties. The proposed technique of simplifying 
nonlinear processes in inhomogeneous media is therefore not restricted to the presented 
example, but covers a wide application range. 

5. Conclusion 

We proposed a theoretical framework for the incorporation of nonlinear effects within the 
concept of transformation optics (TO). In this context, we derived a general expression 
for the calculation of linear and nonlinear electromagnetic material parameters under 
arbitrary coordinate transformations. The transformation law is formulated in a 
manifestly covariant form that allows the simultaneous treatment of electric, magnetic 
and magneto-electric cross-coupling terms and is applicable to both temporal and spatial 
t r ansf ormat ions . 

As a first application example, we calculated the linear and nonlinear constitutive 
material relations in a moving medium. As expected from the transformation of 
electromagnetic fields, the movement of the medium leads to a mixing of the electric and 



(55) 
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magnetic response, where the obtained expressions for the linear constitutive equations 
were in agreement with the relations given in the literature. In addition, we showed that 
the movement of the medium also implies a mixing of the nonlinear material properties 
which describe cross-coupled nonlinear interactions between the E- and £?-fields and 
the medium. This means, for example, that a Pockels medium at rest can display a 
Faraday effect if the material is moved relative to the observer, and vice versa. 

In the final part of the paper we focused on time- independent, spatial coordinate 
transformations which are of particular interest for the design of nonlinear TO devices. 
As an illustrative example of such a device, we presented a twisted nonlinear field 
concentrator and calculated the second harmonic wave that is generated when the 
concentrator is illuminated by a strong laser field. In this respect, we demonstrated 
that sophisticated nonlinear phenomena in complex media can take a much simpler 
form if an appropriate coordinate transformation is applied. 

The considerations have shown that the incorporation of nonlinear susceptibilities 
in the TO approach provides a promising computation method for calculating nonlinear 
effects in moving or inhomogeneous media and offers new opportunities for the design 
of novel optical devices with tailored nonlinear properties. 



Appendix 

Summary of the linear terms in equation (12) 

The linear terms in (12) (the free-space contribution and the linear response of the 
material) can be written as 

V = F + 4iTY aK F 

^ \xv - 1 jiv ~ A.fiu ere 

= \ - C) F ** + ^F aK 

= Q fr K - C) + 47r ^) F m . (57) 
And finally, the last line of (12) follows after the re-definition 47rx^ (new) : = 



Proof of equation (17) 
By using the identity 



1 1 

F — - (F - — F- ) = - (6™™ - o~™ n ) F 
1 *? 2 2 3 3% n 

= \g- 3 9rFmn = -9%B k (58) 
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the components of the electric displacement field are 



D i = V Qi = ATr X ^F ITK 

= 87r X °oiF 0j + ^xTFmn 

= 8n X ° lE 3 - ATrg^xT B 3 (59) 



TT r, mn T > > 0^^. mn A/ f7K I? 



and, accordingly, the components of the magnetic fields are given by 

1 

= -^gr n X°±F 0l - 2*grx*LF» 

= -togTxSLEj + 2ngrg 3 kl X k LBr (60) 
By comparing the expressions of D; L and Hi with 

A = 4e 3 + an., 

H i = (jE j + (vr 1 )iB j , (61) 
we obtain the relations given in (17). 

Number of independent components in x^f 1 "' an/3n 
Due to the symmetry properties 

ai/3i02/32 _ v ai/3iQ2/32 _ v a2/32ai/3i _ v /3iala2/?2 (fi9^ 

X\1V \-V\l ^jXV XfXli v uz v 

(exemplary for n = 2), there are only 6 independent combinations for each upper index 
pair otiPi, e.g. a^i = 01, 02, 03, 12, 13, 23. Furthermore, in the sequence a.\fi\ • ■ ■ a n (3 n 
all n pairs ((Xi/3i) can be permuted as a whole without changing the value of x^l a • 
Consequently, for fixed lower index pair fiu, the number N of independent tensor 
components is equal to the number of combinations of n labels where each label can 
take 6 values (with repetition). This number is given by: 

For instance, N{\) = 6, N(2) = 21, N(2) = 56, etc. Finally, we can also vary the 
lower index pair \xv for which also 6 independent combinations are possible due to the 
above symmetries. Hence, the total number of independent components of the tensor 
X^l l3l "' anl3n is given by 6N(n). 

Intermediate transformations used in equation (45) 
The first intermediate transformation / is given by 

7 

<fi = arctan(y/:r) (64) 
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with the corresponding Jacobian matrix and determinant: 



dr 


dr\ 




/ x y 


dx 


9y 




r r 


d± 


dcf> 




[ — V. SL 


dx 


dyj 




\ r 2 r 2 



*f=\Z S) = ( \ 11' \^\ = \- ( 65 ) 



r 



The inverse intermediate transformation is given by: 

x = r' cos <f> 

y' = r' sin <f> (66) 



The corresponding Jacobian matrix and determinant are: 

A h = lay' dy> = . , > \ A h\ = r (67) 

\w 4v \ sm( P r coscf) J 
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